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ABSTRACT 



This thesis presents an analytical technique for the 
design of compensation for multi-loop servomechanisms. A 
completely general multi-loop system is represented in a 
standard block diagram. For this system a standard deter- 
minant is formed which uniquely represents each block, loop 
and path in the block diagram. Based on this determinant, 
the effect of all possible compensators on the coefficients 
of the system characteristic equation is found. This infor- 
mation is presented in a Location Matrix. A model character- 
istic equation is synthesized from the static and dynamic 
performance specifications. Comparison of the uncompensated 
characteristic equation to the proper model equation shows 
which coefficients must be altered. The Location Matrix 
will show all the compensator -path combinations which will 
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affect the desired change. The method requires no intuition, 
and so can be applied by persons without extensive training 
in the servomechanism field. A step-by-step procedure for 
application of this technique is included as a summary. 
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OBJECT 

The object of this thesis is to present an analytical 
technique for the design of multi-loop servomechanism 
compensation. This technique is set forth under the 
following topics: 

(a) Formation of the Characteristic Equation. 

(b) The Effect of Performance Specifications 
on the Coefficients of the Characteristic 
Equation. 

(c) The Effect of Compensation on the Coef- 
ficients of the Characteristic Equation. 

(d) The Choice of Compensation. 
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Chapter 1 

INTRODUCTION 

and 

METHOD OF APPROACH 

1.1 Introduction . The problem of compensation design for a 
multi-loop servomechanism is not a new one. Indeed, most 
standard servomechanisms or control systems texts give it 
some mention. Most such works treat the problem as an 
extension of the single-ldop problem, and attempt to apply 
single-loop concepts to its solution. 

The penalties of such an approach are numerous and 
serious. The first such penalty arises from the technique 
of getting the system into a single loop configuration. This 
involves a block diagram reduction or some similar procedure. 
As a result of this reduction, the identities of the original 
blocks, and the physical components they represent, are lost. 
Far more serious than this, in a compensation design problem, 
is the fact that the information about the arrangement of the 
components is lost. The analysis of this reduced system would 
show algebraically what kind of compensation was needed. It 
would be impossible to transform this information into a given 
compensator fed around certain main path elements. The whole 
sense of element and path configuration no longer has any 
correlation with the algebra. 

There are two ways out of this dilemma. The designer 
can insert all possible compensators, in all usable paths, 
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and re-do the analysis until he finds one combination that 
gives the desired result. Obviously this is a monumental 
task for any system of reasonable complexity. There are 
graphical techniques which eliminate some of this work, but 
they still involve trial and error. The alternative is the 
intuitive approach. Here the designer must have such vast 
experience that he generally knows the effect of all com- 
pensators, and so can eliminate from consideration many 
which could not work. 

In this thesis an analytical technique is presentee? which 
requires no intuition. It examines the effect of all possible 
compensators without the accompanying drudgery. 

1.2 Method of Approach 

A general multi-loop control system, consisting of all 
possible direct, feedback and feed forward paths, is formulated. 
Its components are presumed to be describable by ordinary, 
linear integro-diff erential equations, with constant coef- 
ficients. When inputs to this system are mentioned, the 
terms "step", "ramp" , and "parabola” of input will be used. 

A step of the input quantity is that input which yields zero 
error when there is unity feedback from output to input. For 
a Type I positional servo, a step would be a step of position. 
For a Type II velocity servo, the step would be a step of 
velocity, etc. By treating such a general system, the ensuing 
results can be applied to nearly all practical systems. 

This analytical compensation design technique uses as its 
cornerstone the characteristic equation of the system. After 
a brief discussion of the methods for choosing direct path 
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components, a method for forming this characteristic 
equation is described. The system is represented in a 
standard block diagram. Using circuit analysis techniques, 
the interaction of the blocks, path and loops in this dia- 
gram is shown in a characteristic determinant of the system. 

It is from this determinant that the characteristic equation 
is written. 

The next step in the procedure relates performance speci- 
fications to the coefficients of the characteristic equation. 
Static specifications such as "following error" and output 
in response to a load torque, place limits on coefficients 
or on ratios of coefficients in the characteristic equation. 
Because of the vast amount of work which has been done on 
simple second order systems, most dynamic specifications 
are expressible in terms of the parameters of such systems. 

The compensated system is made to have a response not sig- 
nificantly different from the specified second order response. 
This is done by making the quadratic characteristic equation 
of the second order system a factor of a model characteristic 
equation for the compensated system. All other roots of this 
model equation are forced to be so large that the second order 
roots dominate the response. 

The effect of all possible compensators on the coeffi- 
cients of the characteristic, equation is found through the 
use of the characteristic determinant. All possible paths in 
the system are uniquely represented in this determinant. 
Closing any of these paths with a compensator adds a cofactor 
to the existing determinant. A rapid method for determining 
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the effect of this path closing on the coefficients of the 
characteristic equation is presented. This information is 
displayed in a Location Matrix, 

At this stage of the design process, the designer has a 
model characteristic equation representing the specifications, 
and an actual characteristic equation representing the uncom- 
pensated system, A comparison of these shows which coefficients 
must be altered by compensation. The Location Matrix indicates 
what compensator, in what location, will effect the required 
change. 

This thesis confines itself to the problem of assuring 
dominance of a desired second order response. The technique 
itself is not so limited, and could be used to force a system 
to perform like any desired third or higher order system. 
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CHAPTER 2 

BASIC DESIGN CONSIDERATIONS 
1. Introduction 

Some servomechanisms perform a computing function. 
Examples of this type of servo are found in fire control 
systems. In these servos, the output or outputs may be com- 
plicated mathematical functions of several inputs. The mathe- 
matical relationships to be instrumented will dictate the 
choice of equipment to be used. The type of servo of primary 
interest in this thesis, however, is one which drives a load 
so that its action duplicates the action represented by the 
input signal. An example of thisL type of servo is the drive 
system for a gun. Its input is a voltage or shaft rotation 
which represents the desired orientation of the gun barrel. 

It is the function of the servo to control the power that 
drives the gun mount so that the desired gun barrel orien- 
tation is achieved. The preliminary choice of equipment to 
perform this type of function is considered in this chapter. 
These considerations result in a basic system with sufficient 
power and accuracy to meet the fundamental requirements set^ 
by the specifications. Subsequent chapters deal with modi- 
fying the basic design to meet specifications concerning the 
more detailed nature of the servo’s response. 
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2. Components of a Typical Servomechanism 

Figure II-l is a diagram of a typical servomechanism. Its 
purpose might be to control the angular position and/or velocity 
of a shaft. The shaft is represented by the block labeled 
'’Load" . The input or command might be the angular position 
of a shaft from a computer. The command receiver and the out- 
put measuring device convert the input and output into quanti- 
ties which can be easily compared. The comparator subtracts 
these quantities and produces a signal proportional to the dis- 
parity between input and output; i.e., error. Since this error 
may be a small quantity, it must be amplified. It is then' 
applied to a controller which actuates the prime mover. The 
prime mover drives the load through the gear train so as to 
null the error. 
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Figure II -1 

A servomechanism can be considered to consist of paths 
through which information flpws. Those paths which conduct 
information in the direction from input to output are termed 
forward paths. Those which conduct in the opposite direction 
are termed reverse paths. In essence, the forward paths of a 
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servo transform the input into the output. The reverse, or 
feedback, paths insure the proper degree of correspondence 
between input and output. 

In a relatively simple system, there is usually one for- 
ward path. However, other forward paths may result from com- 
pensation considerations. To avoid ambiguity, the forward 
path which performs the primary function of changing input tq 
output will be designated the main path. Other forward paths 
will be designated feed-forward paths. 

3 • Some Considerations Governing Choice of Components 

The need for a servo implies that the element which generates 
the input has insufficient power to control the load. Since 
the purpose of the servo is to supply the needed power, choice 
of the prime mover is a logical starting point in selecting 
main path components. The characteristics of the load will 
be known. These characteristics may be expressed in terms of 
its inertia, the viscous friction to which it will be subjected, 
external torques which may act on it and the like. The speci- 
fications will also state the velocities and accelerations re- 
quired of the load. This information will set a lower limit 
on the power and torque capacity required of the prime mover. 

The power supply available will further limit the choice. Other 
requirements on the prime mover will arise from cost, weight 
and size limitations, life expectancy, etc. With these limits 
in mind, the designer must then select a prime mover from those 
commercially available or else design his own. 

Choice of the prime mover will normally place limits within 
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which a gear ratio must be chosen. Ahrendt 1 discusses the 
various factors which influence this choice. Such factors 
as low cost or smoothness of output at low velocities may 
place conflicting demands on the choice of gear ratio. Also, 
the inpedance matching aspect of gear ratio selection may be 
of great importance. In order to choose the optimum gear ratio, 
the designer must have a thorough-going knowledge of the speci- 
fications as well as of the characteristics of the prime mover. 

Choice of the prime mover will probably narrowly limit 
choice of the controller. If a variety of controllers is 
available, the considerations of fineness of control versus 
cost, weight, etc. must be brought to bear. 

Insufficient information will be available in the speci- 
fications to allow choice of the amplifier to be used. The 
gain to be provided will be determined by compensation considera- 
tions. This will not be known until later in the design process. 
When it has been determined, the amplifier may be chosen. In 
addition to the conditions imposed on amplifier selection by 
such factors as available power supply, noise generation, and 
operating environment, the designer must bear in mind the impor- 
tance of linearity. Although sufficient band width is not diffi- 
cult to obtain, saturation effects may present a problem* It 
is uneconomical for a system* s amplifier to saturate before its 
prime mover does, as this means unusable power has been purchased. 

The choice of comparator, measuring device and command 
receiver will be a question of obtaining the best linearity and 
accuracy within the limits set by cost, weight and size speci- 
fications, These devices are described in detail in several 
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12 3 

servomechanism texts. 9 * 

Determination of Component Transfer Functions 

Having chosen the main path components, their transfer 
functions must then be determined. For components which have 
their output independent of the frequency of the input, this 
merely involves determining their sensitivity,, It is desirable 
that this sensitivity will be constant over the range of interest 
in the problem. If sensitivity is not constant, it must be 
assumed so in the remainder of the problem by assigning it a 
representative value. 

The transfer functions of frequency dependent components 
may be determined mathematically by application of basic a-c 
circuit theory and the Laplace Transformation. This technique 
is described in detail by Thaler^. If the electrical parameters 
of the component are unknown, it may be necessary to conduct 
tests on the equipment to determine its transfer function. 

These may take the form of frequency or transient response tests. 
The results of these tests will permit derivation of the trans- 
fer function. Examples of this method are presented in refer- 
ences 5 and 6. 

Having obtained the transfer functions of the equipment 
to be used, it is possible to express the system In this stage 
of its design by a signal flow chart or standard block diagram. 
This phase of design is discussed in the next chapter. 
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CHAPTER 3 

FORMATION OF THE CHARACTERISTIC EQUATION 

3.1 Introduction ,* The characteristic equation is an equation 
in s domain which represents the integro-dif ferential equa- 
tions in time space which describe the flow of information 
7 

through a system* The roots of this equation allow the 
determination of the transient performance and the stability 
of the physical system so described. The purpose of this 
chapter is to develop a standard method of forming this 
characteristic equation for a multiple-loop servo system* 

3*2 Standard Block Diagram * A set of transfer functions are 
developed, as described in Chapter 2, for elements of the 
servomechanism. The block diagram portrays the interconnec- 
tion of these various elements to form the overall system* 

The very nature of the block diagram assumes that the transfer 
function of the tandem conbination of two blocks is the product 

g 

of the individual transfer functions* The overall perfor- 
mance of the system is found by combining these blocks 

9 10 

according to the standard rules for their manipulation . 75 
Unfortunately for the designer, the combining of the blocks 
destroys their individuality and the identity of the mechanism 
or circuit they represent* What is needed is an alternate 
representation of the system which pictures it in more detail 
than the reduced (combined) block diagram* This could take 
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a new form, as in the signal-flow diagrams originated by 
S. J. Mason, 1 '*' or it could be a more ordered manipulation of 
the block diagram, as proposed by Yaohan Chu. 10 It is to 
the latter method that we now devote our attention. 

Certain advantages accrue from standardizing the block 
diagram. In this regard the following symbols are defined: 

a. The Nodal Point or Node : These are summation points 

in the main path of the system. Each node may receive a 
number of signals, but delivers only one. The nodes will be 
lettered in sequence, a, b, c, . . n, from left (input) to 
right (output). 

b. The Pick-off Point : This is a point at which there 

is one incoming signal and two or more outgoing signals. The 
outgoing signals are identical with the incoming one. 

c. The Block Transfer Function : The transfer functions 

in the blocks will be designated by double subscript notation. 
This notation is defined for the three types of blocks as 
follows : 

1) Direct Path Block: These blocks are noted G nN where 

the "n" refers to the node preceding the block. "N" is the 
upper case form of the same letter. 

2) Feedback Block: The feedback blocks are subscripted 

G^ m . The "N” is the same as the upper case subscript of the 
direct path block from which they receive their input. The 
"m" is the designation of the node to which they feed. 



G mN* 



3) Feed-forward Block: Such blocks contain the symbol 

The lower case subscript refers to the main path node 
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preceding the block. The upper case subscript is the same as 
the upper case subscript on the main path block which feeds 
into the same node. 

As examples of the above, consider the following three node 
system: 

Figure 3-1 




Output 



(a), (b), and 0 are nodes. ] — ^ symbolizes a pickoff point. 

Although the system above is represented as having but one 
input and one output, there is no limit to the number of 
either. Internal noise generated in one block may be considered 
as an external input to the node following that block. If 
there is more than one input, each is designated by a Roman 
numeral. If there is a transfer function between an input and 
its entering node it will be subscripted as follows: 




(Input “I”) } 
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3.3 Mathematical Analysis Based on the Standard Block Diagram . 

3.3*1 Formation of the Characteristic Determinant : We 

can obtain a set of simultaneous equations for the system 

represented by a standard block diagram by applying the network 

1*5 

analysis technique of Guillemin. ' This entails writing the 
equation for the signal summation at each node. If we iden- 
tify the signal out of a node with the lower case letter of 
the node, we have the following equations in our example: 

At Node a: 

a = I(G Ia ) - c(G cC G Ca ) - a(G aA G Aa ) - 



At Node b: 

b = a ^ G aA^ “ b ^ G bB G Bb^ “ c ^ G cC G Cb^ 



And, at Node c: 

c = a ^ G aB^ + b ^ G bB^ " °^ G cC G Cc^ 



Rearranging these equations we have: 



(1 + G aA G Aa )a + (G bB G B^ b + (G cC G Ca )c _: (G Ia );C 



(-G 



aA 



(-G 



aB 



)a + (1 + G bB G Bb )t> + < G cC G Cb )c “ 
)a + (-G bfi )b + (1 + G cC G Cc )c = 0 



0 



Inspection of these equations reveals a symmetry of arrarlge- 
ment which would permit writing the equations almost without 
reference to the block diagram,. All the forcing functions are 
on the right. On the left are terms which contain components 
and signals within the system. We can express these equations 
in matrix notation as follows: 
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- 




- 






- 


1 + G aA G Aa 


G bB G Ba 


G cC G Ca 


a 




I(G Ia> 


_G aA 


1 + G bB G Bb 


G cC G Cc 


b 


= 


0 


-G aB 


" G bB 


1 + G cC G Cc 


c 




0 



We call the determinant of the matrix of the coefficients 
the characteristic determinant of the system. To this 
determinant we assign the symbol A * 

3»3«2 Manipulations of the Characteristic Determinant . 
Independent of its physical significance, the determinant of 
the coefficients of a set of linear equations has certain 
properties. These properties are illustrated in the tech- 
nique of equation solving with determinants and in the rules 
for expansion of a determinant. 

The solution of the equations is accomplished as follows: 

1) Replace the column corresponding to the desired 
unknown with the column of forcing functions. The resulting 
determinant, divided by A? is the value of that unknown in 
response to these inputs. 

2) Since the determinant represents a linear system, the 
principle of superposition is valid. Hence, in the case of 
multiple inputs, the output resulting from each may be com- 
puted separately and the individual results added. 

3) The output obtained by replacing the last column of 
the determinant with the input column is the signal from the 
last node. It must be multiplied by the performance function 
following the last node to obtain the system output. 
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The characteristic (homogeneous) equation describing the 
system is obtained by setting the characteristic determinant 
equal to zero. 

The normal rules for expanding a determinant apply to the 
characteristic determinant. In general they require the expan 
sion by co-factors of some row or column. Since we will 
usually be interested in the final output, we should expand 
by elements of the last column. To the co-factor of the ith 
row and jth column we assign the symbol j . Similarly the 
element in the ith row and jth column is called G^. Using 
this notation, the expansion of an "n by n u characteristic 
determinant by elements of the last (nth) column is: 

A = i II i G in A ln 

As an example of the procedure described above for finding 
system output in response to a given input consider the 
following: 



Input to node 11 a 11 is I, Gj a = 1 







-G aA 


1 + G bB G Bb 


Output = 


“W 1 


_G aB 


-G bB 



71 



Output - I (G cC )[G aA G bB + G &b + G aA G bB G Bb ] 

2 

3.3.3 Properties of the Characteristic Determinant . In 
addition to properties inherent in all determinants the charac 
teristic determinant has features peculiar to the system it 
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represents . 

1) The characteristic equation, obtained by equating 
to zero, is the sum of unity plus a number of closed loop 
performance functions. By the nature of the elements which 
comprise the determinant, these functions will be ratios of 
polynomials in s. When this sum is placed over a common 
denominator, its numerator will be the usable form of the 
characteristic equation of the system. 

2) Each element of the characteristic determinant 
corresponds to a path in the standard block diagram. 

3) The diagonal elements of the characteristic deter- 
minant correspond to the paths of inner loops of the nodes. 

The elements above these diagonal elements are all feedback 
paths. Those below the diagonal are feed-forward or direct 
paths . 

4) Each feedback and feed-forward performance function 
appears in only one element of the determinant. Each direct 
path performance function may appear in more than one element. 

We now have at our disposal a standard block diagram. 

This block diagram maintains the integrity of each element 
comprising the system it represents. We form this block 
diagram to aid us in compensating the system. To this end 
we need only include those nodes or pick-offs which represent 
points in the system physically usable for compensation. All 
performance functions between such points should be combined 
into one block to simplify the ensuing algebra. The expansion 
of the standard determinant representing the block diagram gives 
the characteristic equation of the system. Chapter 4 will relate 
performance specifications to this characteristic equation. 
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CHAPTER 4 

THE EFFECT OF PERFORMANCE SPECIFICATIONS ON THE 
COEFFICIENTS OF THE CHARACTERISTIC EQUATION 

1 . Introduction 

The performance specifications state the desired system 
output in response to actuating or disturbing inputs. The 

output is specified in terms of its steady state value and its 

i 

transient behavior. 

The characteristic equation which was evolved in the pre- 
vious chapter describes the system performance in the s domain. 
The standard analysis approach is to find the roots of this 
characteristic equation and transform it to the time domain 
where actual and specified performance can be compared. 

The purpose of this chapter is to relate the* specif ications 
directly to the coefficients of the characteristic equation. 

This eliminates the need for finding the roots, plotting the 
response, and other laborious analysis techniques. 

2. Static Performance Specifications 

The static performance specifications describe the steady 
state value of the system output. In response to an actuating 
input the specification will normally state an allowable error. 
In response to a step we would expect no error. A following 
error will be stated for the response to a ramp input. A 
maximum allowable output is usually stated for the response to 
a disturbing input. These disturbances might be noise or some 
form of load torque. 
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Evaluation of the steady state output to a step, ramp or 
parabola of input position is most useful in studying the effect 
of disturbance signals. As established in Chu’ s dissertation^^ 
the performance function relating the input at any node with 
the output at any other node is 



Output 

Input 



A 



( 4 . 1 ) 



where: i^ node into which input is fed, 

node just preceding output block, 

Gjj ^ performance function of output block 
/\^ j ^ cofactor of the i*'* 1 row and the j*'* 1 
of the standard determinant. 

The output in response to a given input, I(s) is 

1(31 °ijAu 

A 



9 

column 



(i|-.2) 



The steady state value of this output is obtained by applying 
the final value theorem as follows: 



Output = Limit 

ss s -*■ 0 ' 



!(s) OjjAtj' 



A 



C ince /\ = Characteristic Equation 



'A 



( 4 - 3 ) 



( 4 * 4 ) 



where common denominator of the ratios of 

polynomials making up the expansion of /\ . 



Output = Limit 
ss s -> 0 



3 I(s) GjjAij cd A 
Characteristic Eqn. 



( 4 - 5 ) 
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For a step of disturbance, i.e. I(s) = A/s 



Output = Limit 
ss 3^0 



A G 1.tA 






since a^. 



is the limit of the characteristic equation as s goes to zero. 

For a specified maximum output over input ratio, k, where 
k _ output c 



ss 



k> Ll.lt [djjAi] 



0 ^. 6 ) 



A close examination of the equation above will reveal the 

effects of such a specification on the a coefficient of the 

o 

characteristic equation. 

Case I - If k is a number not equal to zero, then a Q is 



specified as being greater than the Limit 

s -*> 0 



g .ijAi,i cd a 



This may or may not limit the size of a Q since this limit may 
be identically zero. 

Case II - If k is specified as zero then the 

Limit G 1T A.. CD a must be zero since a_ cannot be 

s - 0 l J J *3 LX] ° 

infinite • 

In Case I where the Limit j(jj j/\j j CD ^"j is a constant 

we have set one condition on a Q which may be expressed in 
equation form as: 



a^ > N 
o 



Note that 



the Limit G. T /V • CD/\ may be constant containing 
s -+ 0 [_ J J AJ 



variable gains of the uncompensated system. The condition on 
a Q may be met by setting these gains which determine the stiff 
ness of the system. This will then place a lower limit on 
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these gains and restrict their use in future compensation. 

In Case II, if the Limit j CD^J ^ an 

s must be introduced into the numerator of /\^ j which does not 
cancel an s in its denominator; or an s must be placed in 



which is not also placed in the denominator of /\±y This 
will raise the order of the characteristic equation. The former 
may be accomplished by introducing derivative feedback from 
any block before the output block, to the actuating input* s 
node. The introduction of the excess s in is brought 

about by putting integration feedback from other than the out- 
put block to the input node, to which the disturbing inputs 
are applied. 

To illustrate the reasoning behind the statements above 
consider a four node system represented by the following stan- 




The standard determinant for this system is: 



A = 




G 



dD°Da 

0 

0 

1 



For a disturbing input at node "c” we wish the output to 

[ G d D A 3 ^CD A J \ 0. 



be zero; but have postulated that Limit 

s 0 



22 



We need to introduce an s into one of the terms of this express 
ion which will drive the response to zero without affecting the 
response to the actuating input. 

The steady state value of the response to an actuating 
step input, A/s, is (by equation if.5): 



Limit 
s 



l | o ['A ”a] = 



Since and CD^ are common to both expressions a change in 
either will affect both responses. If the s is introduced in 
CD^ it must also appear in the denominator of /\^ so that 
its effect is not felt in the case of the actuating input. It 
must not, however, appear in the denominator of /\->^. The 2nd 
and 4th rows of are contained in both cofactors so they 

cannot be used. The 3rd row appears in /\^ but not in 
The only locations usable for an integration feedback which 
would have the desired result are positions in the 3rd row. 

In our example the only such feedback position is around 

The other method of introducing an s in the expression for 
steady state output in response to a disturbance, without alter 
ing the desired actuated response, is to put an s in the numer- 
ator of /\^ which does not appear in the numerator of /\-^. 
This is accomplished by inserting derivative feedback in a 
first row position. The 1st row is the only one in /V^ which 
is not also in /\^. 

In higher order systems there may be peculiarities of the 



determinant which would allow other positions to be used, or 
which would eliminate some of the above positions. Those 
explained above are the only positions possible if all terms 
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ip the determinant are present. 

Computation of the output quantity is of prime import- 
ance in devising methods of making a system insensitive to 
disturbing inputs. In response to an actuating input we 
might better express steady state performance in terms of 
error since this permits an immediate and direct check on 
the accuracy of the system. If error is defined as the 
difference between the actuating input and its desired out- 
put, then the error equation is: 

E = Input - Output (4.7) 

This may be rewritten as: 

E = 



l _ Output 
Input 



Input. (4.8) 

Expressing Output/input as per equation(4.1) we have: 

^ - G iAi 



E = 



A 



I(s) 



(V.9) 



By the final value theorem, and again expressing /\ as in 
equation (4.4), we have 



Error „ = Limit 
ss s - 0 



sl(s) 



Char.Eqn. 






CD 



Char. Eqn. 



A 



(4.10) 



In response to I(s) = A/s: 
Error 



Limit 
s 



t A f"( 

0 L- 



Char Eqn - G.jA^ 
Char. Eqn. 



i CD / 



] <4 - u) 



Since the steady state error in response to a step will be 
defined as zero we can write: 
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(4.12) 



Limit Ajchar.Fqn. - G ^ jA^ j CD yyjj 

Limit Char. Fqn. 
s -> 0 

Limit Jeharacteristic Equation] = a Q ^ 0, then: 



a ° = cd a1 (4 * 13) 

This defines another condition on a Q , but one which will 
automatically be met if there is unity feedback from the 
output block to the actuating input node. 

In response to an input ramp, i.e. I(s) = B/s : 



Error 

ss 



Limit B 
s -► 0 s 



Char. Eqn. - jAjj CD/\ 
Char. Eqn. 

(4.14) 



For zero error in response to a step, the Limit of the 

s -+ 0 

Characteristic Equation equals the Limit G , T A . , CD a = a . 

s -+ 0 1J ° 

Hence the constant parts of the numerator of equation (4. l4) 
are equal. This means that s is a factor of the numerator. 



Error ,, = Limit B 
s - 0 a 0 



i s + a^ 



g ijA 



cd a (4.i5) 



lAV 



Error _ ai .. 

B ss a“ - ^ . 

o 

This equation puts a limit on the ratio of a-|/a Q in the 
compensated characteristic equation. 

This ratio m^|r not be compatible with the a-j/a Q ratios 
of any of the model equations formed from dynamic specifi- 
cations. The restriction on a-j/a Q may be relieved by makifig 
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the ramp following error identically zero, i.e . , raise 
system to Type II. 

All the conditions on equation order, compensator loca- 
tion and coefficient magnitude which arise from consideration 
of static specifications should be held for simultaneous 
solution with similar conditions yielded by the dynamics of 
the problem. 

C . Dynamic Performance Specifications 

The dynamic performance specifications describe the 
desired transient performance of the system. Since they are 
meaningless in an unstable system, they all imply stability 
even if they do not state it. 

As opposed to the static specifications, where there 
is some uniformity, there are nearly as many ways of express- 
ing a desired transient behavior as there are Process Engin- 

12 

eers. The following summary , while not intended to be all 
inclusive, lists the most common specifications and discusses 
how to calculate them. 

THE TEN MOST COMMON DYNAMIC PERFORMANCE SPECIFICATIONS 



NAME 


TYPE 


MFTHOD OF COMPUTATION 


1 

Gain 

Margin 


Stability specifi- 
cation 


Compute maximum gain 
for stability by Routh- 



Hurwitz criterion. Ratio 
of maximum stable gain 
to actual gain is gain 
margin. 
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2 

Phase 

Margin 



3 

M Peak 



Frequency domain On Nyquist diagram, the 

stability specifi- angle between curve and 

cation. negative real axis at the 

unit circle. On Bode 
plot, find oo c (frequency 
of unity gain), compute 
phase shift at co c , and 
subtract from 180 deg. 

For a well behaved system, 
a phase margin of 45 deg. 
means one overshoot in 
response to a step. 



Frequency domain 
stability spec- 
ification . 



The M curve is the mag- 
nitude of the closed- 
loop response vs. 
frequency. It can be 
calculated from open-loop 
transfer function or from 
values picked off Nyquist 
plot. M peak is defined 
as the peak of closed- 
loop frequency response 
normalized to low fre- 
quency value . 

M = A 

1 - A 

Circles of constant M 
can be drawn on Nyquist 
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4 Stability spec- 

Damping ification. 

Ratio, C 



5 Stability 

Damping Specification 

Factor, K,o* n 
or Decrement 
Factor 



plot. The locus of the 
system must remain out- 
side given M peak circle. 
(An M peak of 1.4 usually 
means one overshoot) 

Can be calculated from 
closed-loop transfer 
function or found directly 
from root-locus plot. 
Damping ratio is defined 
as C in a quadratic term 
of the form: 

s 2 + 2Cco^ + o> n 2 
On the root-locus dia- 
gram, £ equals cos 0 
where © is the angle 
between the horizontal 
and the vector directed 
toward the gain point 
on the locus. 

May be found directly 
from the root-locus. 
Damping factor is defined 
by the factors of a quad- 
ratic system, as was 
damping ratio $ it is the 
real portion of the roots. 
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6 

Percent 

Overshoot 



7 

Bandwidth 



Direct evaluation 
of system 
stability 



Frequency domain 
speed -of -response 
specification 



Where the imaginary por- 
tion of the root deter- 
mines the frequency of 
the damped sinusoidal 
response of a system, 
the real portion deter- 
mines the rate of decay* 

May be calculated from 
response of system to 
step function input* 
Defined as peak of the 
response to a step-func- 
tion input as a percent- 
age of final value* 

Usually defined as the 
frequency at which the 
closed-loop frequency 
response falls to l//2 
its low (zero) frequency 
value* Bandwidth is 
available directly from 
the M curve of the system 
An analogous definition 
would be the crossover 
frequency of the system* 
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8 Speed-of -response One of the simplest 

Rise Time specification definitions of rise time 

is l/bandwidth, another 

is l/( \ peak* Rlse time 
is also defined as time 

to first zero of error* 



9 Overall performance 

"Settle-Out 111 Specification 
Time or 

"Synchronizing"' 

Time 



Settle -out time can be 
calculated analytically 
by the inverse Laplace 
transform, or by picking 
values off the root locus 
plot. This specification 
is most convenient if an 
analog computer is avail- 
able. Settle-out time 
is defined as the time 
required for system 
response to a step-func- 
tion input to approach 
and remain within a 
given tolerance of the 
final value. It is 
usually necessary to 
specify maximum over- 
shoot and steady-state 
error in addition to 



settle -out time 
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10 



ITAE or Integrated 
Value of the 
Product of Time 
and Absolute 
Value 



Overall 

Performance 

Specification 



The ITAE specification 

is defined by the integral 
oo 

J t|e|dt 
0 



Minimize the value of 

this integral for optimum 

system performance. This 

specification yields a 

number that depends on 

system parameters. It 

would probably be an 

improvement to generalize 

it by normalizing to the 

rated output quantity. 

The dimensions would then 
2 

be sec . Another poss- 
ible generalization 
would be to integrate 
to some arbitrary large 
number, such as 10 time 
constants, rather than 
to infinity. This would 
permit the specification 
to be applied to responses for 
which the system yields 
a finite final error. 



ITAE is one of several 
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attempts to define an 
overall figure of merit 
for system operation. 
Rather than simply sum- 
ming the total error of 
a system, the error is 
progressively weighted 
more heavily as time goes 
on. This puts a premium 
on rapid, accurate settle- 
out, and allows for 
unavoidable initial 
large error. 



Due to the complexity and variety of responses of sys- 
tems of order higher than second, very little information 
is available in the literature concerning the characteristics 
of their response. However, the response of a second order 
system is easily computed, fairly simple, and can be charac- 
terized by two parameters. If the system has an oscillatory 
response, the parameters most commonly used are C and oo n , 
the damping ratio and undamped natural frequency respective- 
ly. The use of these two parameters is quite natural since 
the coefficients of the system quadratic equation are simple 

functions of K, and co . 

n 

As seen above, dynamic specifications will not always 
be given in terms of C and cc> n , but translation to this form 
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is possible. As an example, Instrument Engineering . Vol. II ~^ 
provides graphs of calculated second order response in both 
the time and frequency domains which may be used for this 

ILl 

purpose. Chu gives several rules for translating time 

9 

domain specifications into C and co n and Thaler gives others. 
Thaler also includes simplified graphs giving relationships 
between frequency domain specifications and t and oo n as 
well as rules of thumb relating phase and gain margin to 
height of resonant peak. 

Translation of specifications into C and o) n or relatable 
parameters may be other than precise. Analytical relation- 
ships may not exist, as in the case of phase and gain margins. 
Specifications may put varying demands on the physical com- 
ponents of the system such that time dependent coefficients 
are required. They may prescribe a performance unrealizable 
by an order lower than third. In cases such as these, 
reversion to trial-and-error design techniques may be dic- 
tated. If the response of a higher order system can be made 
to closely approximate a second order system, then adequate 
correlation between response and the C, co n type specifications 
is possible. It is to the problem of assuring dominance of 
a given quadratic factor that we now turn our attention. 

A quadratic factor is said to have "dominance” if the 
response to the higher order system closely approximates that 
of the specified quadratic. A high "degree of dominance" 
means a close approximation, i.e. small deviation of actual 
from second order response. The use for which the system 
is intended and the strictness of the specifications will 
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indicate the degree of dominance required. 

In an oscillatory second order system, the correlation 
between coefficients of the characteristic equation and 
the C and co n specifications is 

s 2 + 2C(0S + 00 2 = 0 (4.16) 

n n 

If our system has a cubic characteristic equation with 
one pair of oscillatory roots it will appear as follows on 




where a^, the dominance factor, is defined as shown in 
figure 4-1. The model characteristic equation in terms of 
Bp £ and co n is: 

(s+c^oo^) ( s 2 +2Co) n s+co n 2 ) = 0 

)o) i s 2 *-(2£g,+1)go 2 s+a,co ^ = 0 (4.1 7) 

■in l n in 

The degree to which the response of this system approximates 
a second order response depends on the size of a 1 . A detailed 
examination of the effect of on the deviation of a third 
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order response from that of a second order is shown in 
Appendix A. 

A system with a fourth order characteristic equation 
might take either of the following s-plane configurations: 





Figure 4-2b 



The model characteristic equation in the case shown in 
figure 4-2a is: 

(s+a^CO )(s+a o C0 )(s 2 + 2 CC 0 M S+W 2 ) = 0 
in ti n n n 

s l+ +(2!;+p ;L )(o n s 3 +(l+2CP 1 +P 2 )to n 2 s 2 +(p 1 +2!:p 2 )tt> n 3 s+p 2 (O n l+ = 0 

(4. IS' 



For the case in figure 4-2b the model characteristic 
equation is: 

(s 2 +2a n C 0 co s+a., 2 )(s 2 +2^ 1 co s+oo = 0 

i^nin inn 

s^ ( 2 C C 2 ) a) n s '^ + ( l + 4 ^ ^ 2 + a l 2 ^ °° n 2 s2+ ^ a 1^2 + ^ ^ 1 Q l^ ^ ^n^ S 

+a 1 2 0> n l+ =0 (4 ./<?) 
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In choosing the a's for the case in figure *+-2a a 
reference to the curves of Appendix A is helpful but not 
exact. Each a must be sufficiently large so that its 
effect alone would be within the tolerances allowed. A 
conservative choice of the individual a’s is necessary to 
insure against the possibility of accumulating error of the 
same sign. Determination of the effect of more than one 
real root is a fruitful area for further work. 

In the case where there are two pairs of oscillatory 
roots, there is a slightly different problem. Here it is 
not only the decay of the non-dominant pair, but the size 
of the coefficient of their sinusoidal contribution to the 
time response, which determines their effect. In the 
extreme case, where they are on the imaginary axis, this 
coefficient is approximately l/a times that of the domi- 
nant pair. 

The extreme condition described above also serves to 
place a lower limit on the size of the coefficients of the 
characteristic equation in order to have stability. If 
= 0 the minimum value of the coefficients of the fourth 
order characteristic equation are defined as follows: 

s 4 +(2i:, )eo„ s 3 +(l+a. 2 )cD 2 s 2 +(2C,a, 2 )oo 3 s+a ,V k = 0 (4.2 0) 

in In lin in 

Similar examination of higher order characteristic 
equations shows the relationship between a's, C and o> n , 
and their coefficients. Appendix B lists the various forms 
of these coefficients, including minimum stability values, 
for model equations up to seventh order. 



36 



As a result of the considerations of static and dynamic 
specifications we can now write model characteristic equations 
for a system which will exhibit a specified response* Con- 
siderations of chapters II and III have given us an actual 
characteristic equation which must be forced into this 
specified configuration by compensation* 
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CHAPTER 5 

EFFECT OF COMPENSATION ON THE COEFFICIENTS 
OF THF CHARACTERISTIC EQUATION 

5.1 Introduction . In Chapter 3 we formed the characteristic 
equation which represents the components in our servomechanism. 
The static and dynamic specifications were related to coeffi- 
cients of model characteristic equations, and these are tabu- 
lated in Appendix B. The third phase of our compensation 
procedure consists of determining the effect of compensation 
on the coefficients of the characteristic equation. This effect 
can be considered under two headings: 

a) The Effect of Compensator Location 

b) The Effect of Compensator Type. 

The purpose of this chapter is to show how to quickly deter- 
mine each of these effects. The relationship between them 
and the coefficients will be examined. A notation will be 
introduced which permits rapid evaluation of the effect of 
all possible compensators, without repeated expansions of the 
characteristic determinant. 

The separation of the two effects is possible largely 
because of the unique manner in which the addition of a com- 
pensator alters the characteristic determinant. Before con- 
sidering the two effects separately we will examine the effect 
of compensation on the determinant of a general, four-node 
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uncompensated servomechanism. 



5.2 The Effect of Compensation on the Characteristic Deter- 
minant. In studying the effect of a compensator on the charac- 
teristic determinant, we will use the four node system in 
figure 5*1. 




Figure 5-1 

The characteristic determinant for this system is 
10 0 



A = 



-G 

0 

0 



aA 



bB 



cC 



G dD G Da 



0 

0 

1 



The characteristic determinant is single valued, and so its 
value is independent of the method used in expanding it. 

For example: When expanded by cofactors of the fourth column 

1 0 



-°aA 1 
0 -G 



A= 


“ G aA 


1 


0 


+(-d 8 (i) 




0 


“ G bB 


1 






0 


0 


“ G cC 





bB 
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A = G dD G Da G aA G bB G cC +1 



A similar expansion by cofactors of the third row yields: 



A = (-D 5 (-G bB ) 


1 

- G aA 


0 G dD G Da 
0 0 


+ (-d 6 (i) 


1 

" G aA 


0 G 

dD Da 

1 0 




0 


1 — 1 

o 

o 

0 

1 




0 


0 1 



A = G bB G dD G Da G cC G aA + 1 



Note that both expansions give the same result, and indeed, 
any expansion would give this result* 

If we were to add a compensator to any path previously 
open, we would put a transfer function in a determinant posi- 
tion previously zero. As an example, consider closing the 
path from the output of block G^g to node "a" with an element 
whose transfer function is G g a * The characteristic deter- 
minant of this modified system is: 



An = 



aA 



G bB G Ba 

1 

_G bB 



0 

0 

1 



G dD G Da 

0 

0 



0 




0 


-G cC 


1 






Expanding by cofactors 


of 


the 


first row: 










1 


0 


0 




" G aA 


0 


0 


A m = (-i) 2 d) 

— k m 


- G bB 


1 


0 


+(-l) 3 ( G bB G Ba ) 


0 


1 


0 




0 


_G cC 


1 




0 


~ G cC 


1 



+ • • 



tfO 



(°dD G Da ) 



-G 



aA 
0 -G 



bB 



-G 



0 

1 

cC 



= 1 + + G^G^.G^G^G. 



J bB Ba aA u dD u Da u aA u bB u cC 



This is merely the old A with the term G bB G B a G a ^ added. 

A similar effect is observed is an existing path is 
paralleled. Suppose an element whose performance function 
is G^ was put in parallerl with G aA - This would make the 
transfer function between nodes u a li and ,l b‘ r equal (G qA + GJ^) 
Since no new position in the determinant was filled, the 
expansion remains: 

An = (G aA + G ^A )G bB G cC G dD G Da + 1 

An = G aA G bB G cC G dD G Da + 1 + G aA G bB G cC G dD G Da 

This again is the old A with a new term added. We can gene- 
ralize the above observations as follows: 

A» ' a * <5 - l > 

where G^ is the element added to the i^ h row, column of 
the determinant as a result of adding the compensator. j 

is the signed minor (cofactor) of this determinant position. 
It is the effect of G^A^ j on the coefficients of the charac 
teristic equation that we must determine. 

All the terms in equation (5.1) contain ratios of 
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polynomials in "s'*. These polynomials come from the products 
of terms G N m > etc, in the characteristic determinant. 

These transfer functions represent components of the servo- 
mechanism. By their physical nature, components normally 
contribute a time lag to the system. This lag appears in 
the transfer function as an (s+a) factor in the denominator. 

The expansions of /\ and consist of unity plus a series 
of products of such transfer functions. Because of the pre- 
valence of lagging components these products will have nume- 
rators of order less than (or at best equal to) their denom- 
inators. When the terms comprising /\ are placed over a common 
denominator, the numerator so formed is the usual form of 
the characteristic equation. The common denominator will be 
the product of all the non-repea'ted (s+a) factors, and the 
highest power of s in any of the individual product's denom- 
inators. The numerator contains unity times this common 
denominator. It is this contribution to the numerator that 
contains the highest power of s, and thus sets the order of 
the characteristic equation. 

With these thoughts in mind, consider again equation (5«1) 




(5.D 



Modified Char. Fan. _ Chr~ ” — 



Modified Char. Fan 




(5.2) 




Modified Char. Eqn. = [ C,P, Am] Char. Eqn. + 

C.D.^ 



ij^ij ^ m 
(5-3) 



From this we see two types of changes to the coefficients of 



h2 



the characteristic equation: 

a) If is not unity, i.e., if the compensator 

added new (s+a) factors or a higher power of s to the common 
denominator, the old characteristic equation will be multi- 
plied by the product of these added terms. This will raise 
the order of the characteristic equation. 

b) The coefficients of s k in jCD^ ^ will be 

added to the coefficients of s in CD^m I characteristic 

L CD ^ J 

equation. 

Had the purpose of this chapter been merely to show the 
effect of a compensator on the coefficients of the charac- 
teristic equation, it would be fulfilled. One would only 
have to apply equation (5*3) for each compensator-path 
combination. For an ,l n“ node system and "m w possible com- 
pensators, this requires (n x n x m) applications of equation 
( 5 * 3 ) • This is a process that requires many, many hours of 
algebra. It usually Just isn't done. The designer relies 
on his previous experience to help him pick a usable path and 
is likely to overlook one which would have done the Job 
easier. The remainder of this chapter will develop a tech- 
nique for examining all possible paths in a reasonable time, 
and for displaying this information in a useful manner. 

5*3 The Effect of Compensator Location * We have seen 
(equation 5*1) that adding a compensator has the effect of 
adding a term, jA^ j « to the characteristic determinant. 

Factoring the performance function of the compensator from 
the term G^ A^ j will leave the effect of closing the path 
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in which the compensator was inserted. It is this effect 
that we call the "Effect of Compensator Location". G ir£ir 

G c 

It is the term which would appear as characteristic deter- 
minant element i-j, if the path represented by that element 
were merely closed. 



G ij = G L’ if G c = 1 

Rewriting equation (5*3) for this condition: 



Modified Char. Eqn. = [cD/^* ] Char • E< l n * + G C G lAij CD Z\ m 



(5.4) 



Since there is no denominator in the performance function 
G c = 1, there are no new factors added to the common denomi- 
nator by the inclusion of j . CD^ m = CD^ . 

Modified Char. Eqn. = Char. Eqn. + G^A^ jCD^ (5*5) 

Modified Char. Eqn. - Char. Eqn. = G^A^ jCD^ (5*6) 

(Modified Char. Eqn. - Char. Fqn.) will be a polynomial in s 
containing only those powers of s whose coefficients were 
changed by closing the path. If we could quickly determine 
the powers of s in G^A^ jCDa , for all possible paths, we 
would know what coefficients of the characteristic equation 
were altered by closing those paths. 

To accomplish this we will introduce a new notation. 
This notation will describe the powers of s contained in the 
polynomials which comprise the terms in equation (5*6). 
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Since these polynomials in s arise from the multiplication 
or division of factors (s+a) or ”s” , they will contain all 
the successive powers of s between the highest and lowest 
power present. The symbol so used will be called an "Order 
Number”, and will consist of two numbers, one Arabic and 

one Roman, i.e. (A,R). The highest power of s in the poly- 

A -4-R P 

nomial is s , and the lowest power of s is s . All powers 

of s between these are present. For example: 

4^2 

If (Order of P(s)) = (3,1), then P(s) = a^s +a^s +a 2 s +a^s. 
The algebra of these Order Numbers is much like that of 
exponents. When two polynomials are multiplied, their Order 
Numbers add to give the Order Number of the product. If 
one polynomial in s is divided by another, the difference 
in their Order Numbers is the Order Number of the quotient. 

Order Numbers add and subtract as follows: 

(Ai,Ri> + (A 2 ,R 2 ) = (A lt A 2 , 

As an example of this use of Order Numbers, consider: 

2 ^ 2 
(s+a)(s +bs) = s J + (a+b)s +abs 

Order Number of (s+a) = (1,0) 

2 

Order Number of (s +bs) = (1,1) 

(1,0) + (1,1) = (2,1) 

Order Number of s^ + (a+b)s 2 +abs is (2,1) 

Since the terms of equation (5-6) are polynomials in s of 
the type representable by Order Numbers we can rewrite this 
equation as follows: 

(Coef. Indicator) = (Order of G^) + (Order of j ) + 

(Order of ) 



(5.7) 
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Before examining the use of equation (5*7), it is well 

to become more familiar with these Order Numbers, and to see 

how they are found for the terms of the equation. 

The Arabic number equals the highest power of s arising 

from the combination of the terms of the form (s+a). 

Example Is In (s+a) (s+b) (s+c) , the Arabic number is 3* 

Example 2: In (s+a) (s+b) , the Arabic number is 1. 

(s+c) 

The Roman number equals the net power of s arising from 
terms of the form s 11 . 

Example 3: In s^, the Roman number is III. 

Example 4: In s^, the Roman number is (m-n) . 

s 

As examples of polynomials in which both Arabic and Roman 
numbers occur: 

Example 5: The order of ( s+a ) is (0,1) 

s(s+b) 

p 

Example 6: The order of s (s+a) ( s+b) (s+c ) is (3>H)* 

Example 7: The order of 0 is (- 00 ,- 00 ). To simplify the 

notation, this will be designated (-,-). Note that (A,R)+ 
(-,-) = (-,-). Any sum containing (-,-) has no value. It 
is as if you multiplied by zero. 

The Order Numbers comprising the left side of equation 
(5.7) are obtained as follows: 

(Order of Q^): In determinant positions corresponding 

to feedback paths, is the transfer function of the block 
immediately preceding the feedback pickoff. It would be of 
the form of Example 5 above. In determinant positions corres- 
ponding to feed forward, or parallel to a direct path, it 
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equals unity. The Order Number of a constant is (0,0). 

(Order of A^ j ) : This could be found by expanding the 

minor, and then applying the rules exemplified above. 

A simpler method consists of replacing each existing term 
in the .characteristic determinant with its Order Number. 

The minors are then expanded by the usual determinant rules 
with the following exceptions: 

1) The Order Numbers are added where the terms they 
represent would have been multiplied. 

2) Arabic and Roman numbers are added separately. 

3) The largest Arabic and the largest Roman numbers 
resulting from these additions comprise the Order Number of 
the minor. 

(Order of CDy\ ) : The common denominator of will have 

already been found in the determination of the characteristic 
equation (Chap. 3)* Apply the rules above to this denomi- 
nator as in Example 6. 

With the use of the Order Number technique, a matrix, 
showing the results of equation (5*7) for each path, can 
be formed. This “Location Matrix 11 is constructed as follows: 

a) Replace the existing transfer functions in the Char- 
acteristic determinant by their Order numbers. 

b) Evaluate the (Order of A^j) for each i-J position 
usable for compensation. (Normally all positions) 

c) To this (Order of A^ j ) add (Order of G^) for the 
position concerned. 

d) Add (Order of G^A^) to (Order of CDy\) and display 
the resulting Order Numbers in a matrix. Each Order Number 
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will occupy a position in the matrix corresponding to the 
position of its path in the Characteristic Determinant. 

This matrix will indicate which coefficients can be 
changed by closing a given path. 

To clarify the ideas presented in this section, the 
location matrix for a sample problem will be formed. Con- 
sider the system of Figure 5-2. 




This system has the following characteristic determinant: 
10 0 

s ( s+a ) 

-K, 100 



A 






1 

-K 2 

(s+b) 



Expanding this in terms of elements of the fourth column: 

K-jK^Cs+c) 

A = 1 + s( s+a ) ( s+b'yCs+'d'y 



= s(s+a) (s+b)(s+d)+K 1 K 2 K i+ (s+c) 
s( s+a) (s+b) (s+d) 
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Expanding the numerator and introducing standard symbols for 
coefficients: 



Char.Eqn. 

cd a 



432 

s +a- > s''+a 0 s +a-, s+a A 
•s 2 1 o 

s( s+a) (s+b) ( s+d) 



Applying the rules for formation of a Location Matrix: 



Step a: 



A = 



Step b: A 1X 

^21 

^31 



( 0 , 0 ) 


(-,-) 


(-,-) 


(- 1 ,- 1 ) 


( 0 , 0 ) 


( 0 , 0 ) 


(-,-) 


(-,-) 


(-,-) 


( 0 , 0 ) 


( 0 , 0 ) 


(-,-) 


(-,-) 


(-,-) 


(- 1 , 0 ) 


( 0 , 0 ) 



( 0 , 0 ) A 12 - ( 0 , 0 ) A 13 - ( 0 , 0 ) A^ - (- 1 , 0 ) 
(- 2 ,- 1 ) A 2 2 _ ( 0 , 0 ) A. 23 - ( 0 , 0 ) Z\. 2 lj. - (- 1 , 0 ) 
(-2,-1) A 32 - (-2,-1) A 33 - (0,0) A 3 4 - (-1,0) 



A43. - (-1,-1) A4 2 - (-1,-1) ZV3 - (-1,-1) A41+ - (0,0) 



The Order numbers of the G^'s are as follows: 

G G G G 

L 11 - (0,0) L 12 - (0,0) L 13 - (-1,0) L l4 - (-1,-1) 

G G G G 

L 21 - (0,0) L 22 - (0,0) L 23 - (-1,0) L 24 - (-1,-1) 

G G G G 

l 31 - ( 0 , 0 ) l 32 - ( 0 , 0 ) l 33 - (- 1 , 0 ) l 3 V - (- 1 ,- 1 ) 

G G G G 

^1 - (0,0) \2 - (0,0) L 43 - ( 0,0) V, - (-1,-1) 



Step c: (Order of G^A^ j ) - position in the array denotes 

subscripts. 



(0,0) 


(0,0) 


(-1,0) 


(-2,-1) 


(-2,-1) 


(0,0) 


(- 1 , 0 ) 


(-2,-1) 


(-2,-1) 


(-2,-1) 


(- 1 , 0 ) 


(-2,-1) 


(-1,-1) 


(-1,-1) 


(-1,-1) 


(-1,-1) 



(Order of CDy\ ) equals (3?D* 
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Step d: 


(3,1) 


(3,1) 


(2,1) 


(1,0) 




(1,0) 


(3,1) 


(2,1) 


(1,0) 


Location Matrix = 


(1,0) 


( 1 , 0 ) 


(2,1) 


(1,0) 




(2,0) 


(2,0) 


(2,0) 


(2,0) 



Step b, finding the Order of each minor, merely entails 
scanning the auxiliary determinant of step a to find the 
largest Arabic and the largest Roman sums present. The 
presence of numerous terms in an uncompensated system 

makes this a few minutes work. This is in comparison to 
the hours of algebra required to evaluate ^/\ 

every possible compensator. 

The Location Matrix is the sum of the terms of Equation 
( 5 * 7 ) for all possible compensation paths. The information 
in the Location Matrix, modified by the effect of compensator 
type, will indicate changeable coefficients of the Charac- 
teristic Equation. 

5 * 3 * The Effect of Compensator Type 
5.3*1 Introduction 

Examination of the type of compensator yields informa- 
tion on: 

a. Which coefficients of the Characteristic Equation 
will be changed, and 

b. the extent of that change. 

It is the Order Number of the compensator in conjunction 
with the Location Matrix that indicates which coefficients 
will be changed. The performance function indicates the 
magnitude of the change. 
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Table V-l lists the compensator types to be considered, 
with their Order Number and Performance Function. 

Table V-l Elementary Compensators 



TYPE 


ORDER NUMBER 


PERFORMANCE 

FUNCTION 


Sensitivity 


(0,0) 


K 


Differentiator 


(0,1) 


s 


First Order Lead Filter 


(1,0)' 


(sta) 


Second Order Lead Filter 


(2,0) 


(s 2 +bs+c) 


Integrator 


(0,-1) 


1/s 


First Order Lag Filter 


(-1,0) 


l/(s+a) 


Second Order Lag Filter 


(-2,0) 


l/(s 2 +bs+c) 



5.3*2 Determination of Which Coefficients are Changed 

The Order Number of the compensator, when added to the 
element of the Location Matrix corresponding to its path, 
gives a new Coefficient Indicator. As previously defined, 
the (A,R) pair resulting from this addition indicates which 
coefficients of the Characteristic Equation are modified^ 

i,e, > a A+R’ a A -1+R ’ ’ a R* 

There are three situations in which addition of a com- 
pensator will raise the order of the Characteristic Equation. 
If the order is raised, there will be two changes to the 
coefficients. One change will be to all coefficients due 
to merely raising the order. The change dictated by the 
Coefficient Indicator will be an additional change to the 
coefficients of the higher order equation. 

Before examining the quantitative effects of compensation 
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on the coefficients, the three cases where the order of the 
Characteristic Equation is raised will be considered: 

Case I - This situation arises when integrator com- 
pensation is used. If the addition of a com- 
pensator results in one of the terms of 

l + ZL+i:TrL+Z7rL+ 

r r r 2 r r 3 r 

having more integrations than the overall 
uncompensated closed loop performance func- 
tion, the order of the Characteristic Equation 
will be raised. In the above expression, L r 

signifies loop performance function. IT 

k 

signifies products of loop performance 
functions taken k at a time. That the order 
of the Characteristic Equation is raised in 
this situation is apparent if the Character- 
istic Equation is derived from the system per- 
formance function obtained by application of 
Mason's Rule.^ 



A simple example of this effect may be 
seen by considering the following system: 
Block A Block B 
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Where A = 1 



Kg 

s 



-K, 



A = 1 + K 1 K 2 



Char. Eon. _ s + % K 2 



CD 



'A 

t. . K.K 0 = 0 

Char. Eqn. = s + 12 

For this system the Location Matrix is: 

(0,1) (0,0) 

(0,0) (0,0)^ 

Consider integration feedback from the output of 
Block B to Node b: 



L. M. = 



Block A. 



Block B 




Here, A = 



-K, 



1 + K 2 



A = 



1 + K 2 + K 1 K 2 
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Char. Ean . _ s 2 + K]K?s + Kg 

CD A ~ s 

Char. Eqn. = s 2 + K^KgS + K 2 = 0 

Notice that the order of the Characteristic 
Equation has been raised by this combination 
of compensator and path. If the integration 
had been fed back from the output of Block A 
to Node a, a loop performance function higher 
than first order would not have been generated 
and the order of the equation would not have 
been raised. In this case, 

A = 1 + K x + k 1 k 2 

S S 

Char. Eqn . _ s + (K] ) (l+K?) 

CD^ s 

Char. Eqn. = s + K-^d+Kg) = 0 

Table V-l gives an Order Number of (0,-1) 
for this compensator. Adding this to the 
proper element of the Location Matrix gives s 

Fed from B to b: (Coef. Ind.) = (0,-1) 
Fed from A to a: (Coef. Ind.) = (0,0) 
Since the Characteristic Equation contains no 
negative powers of s, the first indicator (for 
B to b) appears to be contradictory. It is 
this very contradiction that yields the 
information on equation order. The order of 
the Characteristic Equation must be raised 
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by the amount that the "R" of the Coefficient 
Indicator is negative. That this occurred is 
borne out by comparison of the Characteristic 
Equations in the above example before and after 
this compensation. In the case of feedback from 
A to a, the “R 11 of the Coefficient Indicator 
accurately reflects that there is no change in 
order . 

Rule I - The order of the Characteristic Equation is 

raised by the amount that the Roman Numeral of the 
Coefficient Indicator is negative. 

Case II - Another order raising situation arises when 

any of the terms of 

1 + I L + ZJTL + ZTTL + 

r r r 2 r r 3 r 

when placed over a lowest common denominator, has 
a numerator of higher order than the uncompensated 
equation. Again, this is apparent if the Charac- 
teristic Equation is derived from the system per- 
formance function obtained by application of 
Mason's Rule.^ 

For example, if the sample system is compen- 
sated as follows: 

Block A Block B 
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then, A = 1 



s 



-K^ 1+ K 2 (s 2 +bs+c) 



s 



A = 1 + KjK2 + K2(s 2 +bs+e) 



s 



s 



Char. Earn _ K?s 2 + (K2b+l)s + (Kpc + KjK2 




s 



+ (K-j+c) = 0 



Adding the Order Number of this compensator, (2,0) 
to the proper Location Matrix element gives a 
Coefficient Indicator of (2,0). Again there is an 
apparent contradiction in that the highest coeffi- 
cient indicated to be changed is not present in 
the uncompensated Characteristic Equation. This 
contradiction indicates that the order of the Char 
acteristic Equation must be raised. 

Putting this compensator in the path from the 
output of Block A to Node a yields: 



The Coefficient Indicator in this case is (2,1), 
which again indicates a raise in the order of the 
Characteristic Equation, the coefficients to be 
changed being a^j an( ^ a l* 



^ = Char. F,an. _ K 1 s 3 + Kjbs 2 + (K-|C+l)s + K-1K2 




Char. Eqn. = I^s 3 + K-jbs 2 + (K-^+Ds + KjKg = 0 
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Rule II - The order of the Characteristic Equation is 

raised by the amount 

(A+R) Coef. Ind. ' (A+R) Order of CD^ 

Note: If an occasion arises where both Rules I 

and II apply, first raise the Order of by 

the amount indicated in Rule I, put "R" of the 
Coefficient Indicator to zero, then, if necessary, 
apply Rule II. 

Case III - A lag filter will always introduce its denom- 

inator in CD^ . This will raise the order of the 
common denominator by an amount equal to the order 
of the filter. Since in practical servos, the 
order of the Characteristic Equation equals the 
order of the common denominator of it too 
will be raised by this amount. 

Rule III - The order of the Characteristic Equation is 

raised by the amount the Arabic Numeral of (Order 
of the Compensator) is negative. 

5 * 3*3 Determination of Quantitative Changes to Coefficients 
The Modified Characteristic Equation, Equation ( 5 » 3 )> 
consists of two terms: 

Mod. Char. Eqn. = CD /\m (Char. Eqn.) + j CD /\m 

This indicates that there are two contributions to the change 
in the coefficients. The first term's effect is predictable 
in the case of simple compensators. It is the amount all 
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coefficients are changed when the order of the Characteristic 
Equation is raised by applying Rules I and III above. 

In Case I, CD/\ m S |R| where the "R“ 

CD A is the Roman 

Numeral of the 

Coefficient 

Indicator. 

In Case III, CP Am = (s+a) or (s 2 +bs+c) 

cd A ~ 

^ depending on the 

order of the fil- 
ter used. 

In these cases, the initial effect on coefficients can be 
readily calculated. 

The change to coefficients a A+R , a A-l+R’ etc *> as indi- 
cated by the Coefficient Indicator can only be found by eval- 
uating -j CD /\ m for the path concerned. Evaluating the 

change on coefficients caused by raising the order of the 
Characteristic Equation in accordance with Rule II can only 
be determined by evaluating 
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CHAPTFR 6 

SELECTION OF COMPENSATOR 
1 . Introduction 

The problem to be considered in this chapter is that of 
choosing the combination of compensator and location that can 
most easily change the characteristic equation of the system 
to one established by the specifications. This requires a 
comparison of the coefficients of the uncompensated charac- 
teristic equation to those of a minimum model equation. The 
actual coefficients which exceed a minimum allowable value 
will be used to select a working model equation. The work- 
ing model equation permits calculation of a compensated char- 
acteristic equation consistent with the specifications. 

There are four equations of interest in this chapter. 

The uncompensated characteristic equation describes the 
servomechanism before compensation. Its coefficients are 
set by the basic hardware. It is these coefficients which 
must be changed. The minimum model equation is the charac- 
teristic equation of a system with a transient response which 
has: 

a. A specified dominant damped oscillatory mode, 

b. A minimum number of first order modes, (none or one), 

c. Non-dominant undamped oscillatory modes. 

No first order modes will exist for even ordered systems. No 
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non-dominant oscillatory mode will exist for a third order 
system. Its coefficients are those of the characteristic 
equation of a system verging on instability. The working 
model equation is the characteristic equation of a system 
with a transient response which has a specified dominant 
damped oscillatory mode. Its non-dominant modes will be 
damped oscillatory and/or first order. Its coefficients are 
those of a characteristic equation of a stable system. The 
compensated characteristic equation describes the servo- 
mechanism after compensation. Depending on its stage of 
development, its coefficients may be purely numerical or 
may contain unspecified parameters. Compensation which 
changes the coefficients of the uncompensated equations to 
those of the compensated equation will result in a system 
with the specified response. 

2 . Choosing the Minimum Model Equation 

There are three factors which determine the choice of 
a minimum model equation. They are: 

1. It must be of the same order as the uncompensated 
characteristic equation. 

2. The quadratic made up of the specified C and co n 
must be a factor of the model. 

3. The dominance factors, a's, must be sufficiently 
large to insure that the response of the model has the 
desired correspondence with the response of the quadratic. 
The minimum model equations in Appendix B are so con- 
structed that substitution of the desired C, co n and a’s 
will give them the properties listed above. 
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3 . Determination of the Compensated Characteristic Equation 

The compensated characteristic equation must have coef- 
ficients at least as large as those of the minimum model. 

If its coefficients are less than those of the minijnum model, 
the system must be unstable. A comparison of the uncompen- 
sated characteristic equation to the minimum model equation 
will show which coefficients must be changed. There are 
two possible conditions arising from this comparison. Either 
all of the coefficients are below the minimum, or a number 
of them are already sufficiently large. In the first case 
substitution of non-dominant damping ratios and/or acceptable 
dominance factors into the working models of the proper order, 
as given in Appendix B, will set the compensated characteristic 
equation. The ratio' of a^ to a Q in the compensated charac- 
teristic equation may be set by static specifications. This 
may limit the choice of a's and non-dominant C's. 

In the case where some or all of the coefficients are 
already above the minimum, several factors must be considered. 
These are; 

1. The interdependence of the coefficients of the model 
equation. 

2. The position of the satisfactory coefficients. 

3. The effect of static specifications. 

The interdependence of coefficients is expressed in the 
working model equations of Appendix B. In an n^* 1 order model 
the n coefficients are functions of the specified £ and o) n , 
and of (n-2) dominance factors and/or non-dominant C's. 

This gives rise to n equations in n unknowns. Two of these 
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unknowns, C and co n of the dominant quadratic, are set by- 
specifications. The remaining (n-2) degrees of freedom 
may be constrained by fixing either coefficients or a's. 

If any (n-2) coefficients of the uncompensated characteristic 
equation are above the minimum values, they will determine 
the required value of the remaining two coefficients. If 
fewer than (n-2) coefficients are satisfactory, the 
remaining degrees of freedom may be fixed by an arbitrary 
choice of a*s. Required values for the remaining coef- 
ficients can then be calculated. These satisfactory and 
calculated coefficients comprise the compensated character- 
istic equation which will require minimum change in coeffi- 
cients. The satisfactory coefficients should first be com- 
pared to the model which contains all real non-dominant roots. 
If they are not sufficiently large to fit this model, then 
the oscillatory working model equation must be used. 

The second consideration arises from the fact that all 
compensators affect successive coefficients in the character- 
istic equation. As developed in Chapter the independent 
coefficient changes are found by evaluating G ij A ij CD Am‘ 

This will be in the form of Ks n (s+a) ( s+b) etc. When 
expanded, it will contain successive powers of s. This fact 
is also revealed by the Coefficient Indicator. Because of 
this fact, the satisfactory coefficients used to determine 
the compensated characteristic equation should be chosen so 
that the coefficients remaining to be changed are successive. 

The static specifications placed lower limits on some 
coefficients of the compensated characteristic equation. 
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These limits may force some coefficients to be larger than 
the dominance factors dictate. The coefficients so affected 
must be included in the (n-2) degree-of -freedom considerations 
above . 

Comparison of the desired Coefficient Indicator and the 
Location Matrix will determine the type of compensator required 
to change the necessary coefficients. If it raises the order 
of the equation, the abpve procedure must be repeated start- 
ing with the minimum model equation of the higher order. The 
uncompensated equation to be compared with the model can be 
formed knowing the order raising effects described in 
Chapter 5* 

*+. Evaluating the Parameters of the Compensator 

The of the required compensator is added to the 

old characteristic determinant. This results in a charac- 
teristic equation in which the coefficients requiring modifi- 
cation contain the parameters of the compensator. Equating 
these coefficients to those of the compensated characteristic 
equation will determine the value of the parameters. The 
compensator must have sufficient degrees of freedom to 
match the independent coefficient changes required. 
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CHAPTER 7 

COMPENSATION DESIGN PROCEDURE 
7*1 Introduction 

The advantages of the technique described in the previous 
chapters are three- First, it involves no graphic presenta- 
tions as required by other methods- Second, it does not depend 
on the intuition of the designer for its choice of compensation 
to achieve the desired result. Finally, it requires little 
of the extensive algebra inherent in other analytical methods. 
Since there is no intuition or experimental graphing involved, 
an ordered analytical approach to any compensation problem 
is possible. 

This method divides naturally into four phases, which 
are the subjects of Chapters 3 through 6. These chapter 
divisions were made primarily in order to explain the tech- 
nique. They also follow in the proper logical order for 
application to the solution of a problem. The step-by-step 
procedure in this chapter thus serves as a summary of the 
paper as well as a guide for the compensation- designer . 

7*2 Step-by-step Guide to Compensation 

Given: I. A servomechanism for which the direct path 

components have been selected. This servo may have 
been partially compensated or may just have unity 
feedback from output to input. 
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II. Static and Dynamic Performance Specifications 
to be met by the system. 

Problem: To design compensation to meet the specifications. 

Solution: 

Phase I: Forming the characteristic equation. 

A. Arrange the system in a standard block diagram. 

COMMENT: Combine the components between usable 

pickoffs and nodes into single blocks. 

B. Form the characteristic determinant. 

C. Write the characteristic equation by setting 
the expanded characteristic determinant equal to 
zero. 

COMMENT: Note the Order Number of CD^. 

Phase II: Relating Specifications to Coefficients of 

Model Characteristic Equations. 

A. Express static specifications in terms of 
restrictions on Sq and a-^ of the compensated 
characteristic equation. 

B. Interpret dynamic specifications in terms of 
second order response characteristics $ i.e., C and 

v 

C. Choose minimum acceptable dominance factors 
(q's and non-dominant C's) to insure second order 
dominance . 

Phase III: Determining the Effect of Possible Compen- 

sators on Coefficients of the Uncompensated Character- 
istic Fquation§ i.e., Forming the Location Matrix. 
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COMMENT: Chapter 5 shows the details of this phase. 

Phase IV: Choosing the Compensator. 

A. Form the minimum model equation using a’s, C’s and 

co found in Phase II. 
n 

COMMENT: All coefficients of the compensated equation 

must be greater than those of the minimum model to 
insure stability. 

B. Choose the Working Model Equation. 

1. Insure compliance with dynamic specifications. 

Put sufficiently large a*s (Phase II-C) and the required 
C and oo n (Phase II-B) in the working model equations of 
the proper order. 

2. Insure compliance with static specifications. 
Compare a^/ a Q ratio, of the various forms of the working 
model equation. Choose for a working model that form 
which, with allowable non-dominant zetas, (Phase II-C) 

has the smallest coefficients that meet the specifications. 

3. Remove the effect of static specifications from 
working model equation. In the event that the minimum 
model equation will not meet static specifications, the 
basic system configuration must be changed by addition 
of a component which reduces following error to zero. 

This will invariably mean the addition of an integration. 
Compensation of this type will introduce a zero in the 
system performance function. Set this zero to a value 
which will remove its effect on the dominant mode of the 
response (See Figure A-2). The design procedure must 

be started again with the revised system. 
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C. 



D. 



Form the compensated Characteristic equation, 

1. Compare the uncompensated characteristic equa- 
tion to the working model, 

2. If n-m-2 coefficients of the uncompensated equa- 
tion are above the working model coefficients and 
leave successive coefficients to be altered, use 
them to fill some of the n-2 degrees of freedom. 

Set m arbitrary a 1 s or non-dominant C*s. Compute 
the m+2 remaining coefficients of the compensated 
characteristic equation. If no coefficients are 
above those of the working model, it is the com- 
pensated characteristic equation. 

Select the Compensator. 

1. Form the desired Coefficient Indicator from the 
coefficients in C-2 above which must be changed. 

2. Compare this Coefficient Indicator to the Loca- 
tion Matrix to find the simplest possible compen- 
sator. The compensator selected must introduce 
enough additional variable parameters so that the 
total number equals the number of coefficients to 
be changed plus the number of zeros of the system 
to be set (See COMMFNT) 



3. Form ^In^l J^-1.1 , which is the system perfor- 
mance function relating the actuating input at 



node "n" to the putput of block G 






This 



expression will contain the adjustable parameters 

•CT 

of the compensator. The compensated characteristic 
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equation will be the denominator of this function. 

4. Equate the parametric form of the coefficients of 

the denominator of the system function to the coef- 
ficients of the compensated characteristic equation. 
Also set the parameters of any numerator zeros to 
a value which gives a satisfactory zero (i.e., 
sufficiently large a). 

COMMENT: If the only possible compensator raises the 

order of the characteristic equation, apply the 
order raising rules of Chapter 5 to the uncompen- 
sated characteristic equation. This order raising 
will put a zero in the system function. Repeat 
phase III using higher order model. 
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CHAPTER 8 
SAMPLE PROBLEM 

8.1 Introduction 

The problem to be worked in this chapter is a modification 

of an analysis problem presented by Chestnut and Mayer The 

problem, as presented in Chapter 13 of the reference, consisted 

of determining the performance of a position control system 

subjected to a load disturbance. Compensation had already 

been accomplished. It consisted of feedback of the first and 

second derivatives of the output quantity. The signal fed 

2 

back was modified bv a filter of the form Ks (s+a) • 

(s+bKs+c7fb+dr 

The specifications were estimated from the performance of the 
system as compensated in the reference. 

8.2 Statement of the Problem 

Given: X The system to be compensated has the primary require- 
ment of position control. However, it must respond properly, 
within specified limits, when subjected to a designated 
disturbance load torque. The main path components consist 
of the following: 
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Component 
Preamplifier 
Amplifier, Amplidyne 

Motor and Load 



Transfer Function 



K 1 not specified 



1500 

(s+31.$) 



850 

(s+17.^) 



4.78 

s(s+9.l) 



Disturbance Load/Voltage 

Transducer K T = .855 volt/ft. lb. 



II Static Specifications 

Steady State Error in response to step 
Following Error in response to ramp 



= 0-0 

= lO^rad/rad/sec 



Maximum output in response to load torque = 10 S'ad/ft-lb- 
Dynamic Specifications: 

Peak Overshoot = 1«2 



Settling Time to 1°/ Q of steady state 



in response to step 
Problem: Design compensation such that the system will 

conform to specifications- 
8-3 Solution of the Problem 

Phase I: Forming the characteristic equation* 



= -6 sec 
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B. 



A = 



1 

-1500 
s+31 . 5 
0 






4.78 

s(s+9.l) 



0 . 

A _ , . (6.09 x 10 ) Ki 

^ 1 s(s+31«5 ) (s+ 17-6) (s+9.1) 

A - s(s+31. 5) ( s+17.6) (s+9.1 )+(6.09 x 10 6 ) Ki 
^ " s( s+31.?) (s+17.6) (s+9.1) 

A _ s tf +58.2s 3 +1001s 2 +5040s+(6.09 x IQ 6 ) Ki 
^ " s( s+31. 5) (s+17.6) (s+9.1) 

Char. Eqn: + 58.2s 3 + 1001s 2 + 5040s + (6.09 x lO 6 )^ = 0 

(Order of CD^) = (3,1) 

Phase II: Relating Specifications to Coefficients of Model 

Characteristic Equation* 

A* Static Specifications 

1* Steady State Position Error in response to step = 0* 

a 0 = f G dDAi4 CD A] 



a^ = Limit 

° s 



t f 4.78 (1.275 x 10 6 )Ki 

0[_s( s+9.1) (s+31. 5) (s+17.6) 



s(s+9.1)( s+31. 5) ( s +17.6) 



a 0 = (6.09 x 10°)^ 



H 



The value of a Q is independent of this specification 
because of the unity feedback* 

—L). 

2. Following Error in response to ramp = 10 sec. 
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3. 



_K 

Error = =10 > 

ss 



5040 



(6.09 x 10 )K X 



K. > S°!+0 — p = 8.28 or f 1 < 10 -4 

1 6.09 x 10 o 

Response to Load Torque = 10 -l+ rad/ft. lb. 




10 "^ > Limit [" (.855)(4.78)(s)( s+17.6) (s+31. 6) ( s+9.1 ) 
s "" 0 [ (6.09 x 10 6 )(s)(s+9.1)K 1 . 

10 _lf > (.85 c O( 1 4-.78)(17.6)(31.5) 

(6.09 x lO 6 )^ 

K x > 3.^3 

B. Dynamic Specifications. 

From reference 13? Peak Overshoot of 1.2 corresponds 
to (.6 > C > •-?). We will choose the conservative 



C = .6. 

From reference 10, Error is reduced to within l°/ 0 
of final steady state value in four time constants, and 



Settling Time = 4 = .6 sec. 

Cco 
^ n 

L. . 

a> n = = 11.1 rad/sec 

C. Choosing Dominance Factors. 

From Appendix A, > 3 and > 2 are acceptable 
values of dominance factors. If non-dominant roots are 
complex, use a > 3 and > 3° 

Phase Ills Determining the Effect of Possible Compensators. 
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Auxiliary /\ = 



Order Numbers 
of -s 



Order Numbers 
of G l 's 



(0,0) 


(-,-) 


(-,-) 


(0,0) 


(-,-) 


(-1,0) 


(-,-) 


(-,-) 


(0,0) 


(0,0) 


(-3,-D 


(0,0) 


(-2,-1) 


(-2,-1) 


(-1,-1) 


(-1,-1) 


(0,0) 


(-1,0) 


(0,0) 


(-1,0) 


(0,0) 


(0,0) 


(0,0) 


(0,0) 



(-,-) 


(- 1 ,- 1 ) 


(-,-) 


(-,-) 


( 0 , 0 ) 


(-,-) 


(- 1 , 0 ) 


( 0 , 0 ) 


(- 1 , 0 ) 


(- 2 , 0 ) 


(- 1 , 0 ) 


(- 2 , 0 ) 


( 0 , 0 ) 


(- 1 , 0 ) 


(- 2 ,- 1 ) 


( 0 , 0 ) 


(- 1 , 0 ) 


(- 1 , 1 ) 


(- 1 , 0 ) 


(- 1 ,- 1 ) 


(- 1 , 0 ) 


(- 1 ,- 1 ) 


( 0 , 0 ) 


(- 1 ,- 1 ) 



(Order of CD^) = (3,1) 



Location Matrix= 



(3,1) 


(2,1) 


(1,1) 


(0,0) 


(0,0) 


(2,1) 


(1,1) 


(0,0) 


(1,0) 


(1,0) 


(2,1) 


(1,0) 


(2,0) 


(2,0) 


(1,0) 


(2,0) 


g the Compensator o 






co = 

n 


For o-,=3 
13.33 






iw 


= 1,233 








= 1^,080 







A 
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a o = a l 2co n 4 = 1 37,000 

s 4 + 13.33s 3 + 1233s 2 + 14080s + 137000 = 0 

Note that a 2 and a 1 of the uncompensated system are 
below minimum values. 

B. a^ for minimum model = #103 
a o 

aq for all complex non-dominant roots = .126 

% 

a^ for all real non-dominant roots = .156 
a 

o 

Since a 1 /a Q from the model equations can never be 
less than .103, the condition on a^/a 0 imposed by static 
specifications must be relieved by making the following 
error independent of this ratio. This effect can be 
achieved by compensating so as to have zero following 
error. This implies making the servo a velocity control 
system. However, it may not lose its position control 
ability. Paralleling the main path with another forward 
path containing an extra integration will result in both 
position and velocity control. This effect can also be 
obtained by placing the extra integration in a feed for- 
ward path from the first to second nodes thereby utilizing 
succeeding integrations in both forward paths simultaneous- 
ly. Using this compensation, the transfer function, 

becomes K^(s + K^/Kl) 6 

“s 

In addition to raising the order of the characteristic 
equation, a zero of the system performance function will 
be introduced at s = -K 2 /K^. Set this zero such that its 



7 ^ 



a > 2, s&y K 2 /K 1 = 25 « It is now necessary to recommence 
the solution. 

Phase I: Forming the characteristic equation. 




B. 



A = 



1 

0 



0 

1 

_ -lSfifi- 

s+31.5 



o 

0 

1 

850 

s+17.6 



4.78 

s(s+9-l) 

0 

0 



^ + 58. 2s 1 * + 1001s 3 + 5040s 2 + 6.09 x 10 6 Ki s + 6.09 x 10 6 Kp 



s (s+9.1 )( s+ 17.6) (s+31. 5) 



Char.Eqn. = s 5 +58.2s l+ +1001s 3 +5040s 2 +6 .09 x 10 6 K 1 s+6.09 x 10 6 K 2 =0 
(Order of CD^) = (3,11) 



Phase II: Relating Specifications to Coefficients of Model 

Characteristic Equations. 

A. Statis Specifications: 

1. Steady State Position Error in response to a 
step = 0 
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2 . 



3. 



a Q = Limit jG dD A llf CD^j = 6.09 x 10 6 K 2 

As before, the value of a Q is independent 
of this specification. 

Following Error in response to a ramp = 0 



Error = Limit 1 
s -+ 0 s 



1 - G. A 






= Limit S 4 + 58.2s 3 + 1001s 2 + 5040s 
s -+ 0 Char. Eqn. 

Therefore, the values of the coefficients are 

independent of this specification. 

-4 / 

Response to Load Torque < 10 rad/ft.lb. 



10“^ > Limit 



10" 4 > Limit (.855X4.78) (s) (s+17. 6) (s+31.5) 
s -> 0 Char. Eqn. 

-4 

10 * > 0 

Therefore, the specification on load torque 
disturbance is met. 



u IId 

Char 



A-4 cd /\1 

.Eqn. J 



B. Dynamic Specifications: 

Same C and o> n are required. 

C. The same dominance factors are acceptable. 

Phase Ills Determining the Effect of Possible Compensators. 

The (Order of G^) for the G^ a feedback position is changed 
'to (1,-1) by the addition of the parallel integration. The 



(Order of CD^ ) is changed 


to (3,11). 








"(2,1) 


(2,11) 


(1,11) 


(0,1) 


New Location 


(0,1) 


(2,11) 


(1,11) 


(0,1) 


Matrix 


(1,1) 


(1,1) 


(2,11) 


(1,1) 




(2,1) 


(2,1) 


(1,1) 


(2,1) 



76 



Phase IV: Choosing the Compensator* 

A. Using a 1 = 3 and a 2 = 2, the minimum fifth order 
model is: 

s 5 + 35.53s 4 + 1527s 3 + 42,100s 2 + 464,000s + 3,030,000=0 
Since this is now a Type II system, it has zero 
error in response to inputs of position and velocity* 

The a 1 /a Q ratio represents velocity error in response 
to an acceleration input* 

The requirement that the working model meet static 
specifications is removed* 

B. 1. Since no working model had to he formed to con- 
sider the static specifications, the initial comparison 
will he made to the minimum model. With any easily 
realizable values for K 1 and K 2 , the a^, a-j_ and a Q of 
the uncompensated equation are above those of the minimum 
model. Since a^ is the only fixed coefficient above 

the minimum model, it will be used as a satisfactory 
coefficient* and K 2 will be left unspecified* 

2. n-m-2=l and m=5-2-l=2 arbitrary dominance 

factors* Solving the parametric form of the a^ 
coefficient of the oscillatory working model, 
with a 1 =3 and a 2 =2, yields? 

(2C-L + 2a 1 C 2 + a 2 )o) n = 58*2 

( 1*2 + 6 C 2 + 2 ) 11*1 = 58*2 

C 2 = o53^ 

This non-dominant damping ratio is larger than needed, 
so C 2 = o3 can be used to resolve for a larger a 2 * 
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(1.2 + 1.8 + a 2 )ll.l = 58.2 
a 2 = 2.24 

Using these values for o's and C's, the remain- 
ing coefficients can be found* 
a 3 = 2820 
a 2 = 54600 
a x = 564000 
a Q = 3415000 

The compensated characteristic equation iss 

s 5 + 58*2s If + 2820s 3 + 54600s 2 + 564000s + 3415000 = 0 
C. 1. The coefficients to be changed are a^, a 2 , 
and a Q . a^ and a Q could be changed independently of 
additional compensation by adjusting and K 2 , but 
K 2 /Ki has been set at 25» Therefore it is necessary 
that compensation provide for varying a^ and a Q indepen- 
dently of the other* In order to keep the coefficients 
to be changed successive, K 2 will be determined* From 
the characteristic equation developed in Phase I-C and 
the compensated equation developed in Phase IV-C-2s 

6.09 x 10 6 K 2 = 3,415,000 

K 2 = .56 

Therefore, = K 2 /25 = .0224 

The 11 semi-compensated" characteristic equation is nows 

s 5 + 58.2s L *' + 1001s 3 + 5040s 2 + 137,000s + 3,415,000=0 
The compensated characteristic equation iss 

s 5 + 58.2s 1+ + 2820s 3 + 54,600s 2 + 564,000s + 3,415,000=0 



7 8 



Comparison of these two equations gives the desired 
Coefficient Indicator to be (2,1) o 

2. The compensator should be the simplest one 
which has sufficient adjustable parameters to change 
the required number of coefficients. The compensator 
must have three adjustable parameters since there are 
three coefficients to be raised. 

The compensator must be in a feedback path since 
feed forward compensators reduce the size of the coef- 
ficients they effect. Comparison of the desired Coef- 
ficient Indicator and the Location Matrix indicates 
that feedback position G Dd requires the simplest com- 
pensator with necessary parameters, one with order 
number (2,0). This can be achieved with a transfer 
function of the form K(s+a)(s+b) which will arise 
from the followings 



to b^Of- 



K-^s 

"v" 



■ from G 



dD 



PF 



k 3 (s +K ‘ +// K 3 s+K ^ // K 3 ) 



An A + G ijAu 



= s ? +‘T8.25Aioois 3 +‘?040s 2 +137.000s+3.4l5.000 
s^(s+9-1) (s+17.6) (s+31«5) 



. (4,78) (1500) <850)K^(s 2 + K 4/K^3+ K 5/K3 ) 
sTs+9 <fl ) ( s+17 « 6 ) ( s+31 • 5 ) 
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VW4-. ■ (1) ( J&s) 



Ia^dD^-14 



G ^a G , 

71 



137.000(3+2$) 
4, ,, , .6., 



s 5+ 58 . 2 s 4 + ( 6 . 09x 10 6 K 3 +1001 ) s 3 



c: 



( 6 . 09x10 6 K 1+ + 5040 ) S 2 + ( 6 . 09x10 6 K ? ) s+3 , 41 9, 000 



Equating coefficients of the two forms of the compen- 
sated equations give the following relationships; 

a^ = 6.09 x 10 6 + 1001 = 2820 or = .000299 

a 2 = 6.09 x 10 6 + 5040 = 54,600 or = .00813 

.6 



a ± = 6 o09 x 10 = 564,000 



or = o0925 



8o4 Checking the Solution 

The compensated system has the following block diagrams 
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An 



1 0 

-(.0224+. 56) 1 

s 



0 



>.78 

s(s+9ol) 



o [AfeJ <• 



299x10 



,-3.2 



+.813x10" 2 s+.0925) 



0 

0 



_ _1500_ 
s+31.5 



1 



0 



8 50 

s+17.6 



0 

1 



^ _ s ? +58.2s lf +2820s 3 +54.600s 2 +564.000s + 3.415.000 

m s 2 (s+9«l) (s+31-5) (s+17. 6) 



For an input of a unit step to node a and output from 
block the system performance function is? 

137,000(s+25) 

PF = s(s+24„9)(s 2 +20s+1109)(s 2 +13.3s+123.2) 

The response of the system in time domain can be found 
by the method of evaluating residues to be? 

e out (t) = l+.0011e“ 2l+ “ 9t +.137e" 10t cos(31.8t-30 o .6) 

+1 .48e ” 6 “ 66t cos ( 8 . 88t-139° .5) 

It is apparent that the contribution of the real root is 
negligible, being only on the order of l°/ 0 at its greatest « 

The contribution of the non-dominant oscillatory roots, though 
somewhat larger, is still less than one-tenth of the contri- 
bution of the dominant roots • It decays quickly and by t = »2 
adds only about 2°/ q to the value of the output <, At t = *6 sec, 
the value of the output is about <, 98 , a 2°/ o deviation from 
final value • A 1°/ deviation was specif iedo At t = <,27, 
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peak overshoot occurs and equals 1*24 5 compared to the sped- 
fled 1 * 20 . Failure to meet specifications exactly in these 
two instances is attributable to placing the zero too close 
to the origin* This indicates a need for further study of 
the effect of zeros upon the transient response, or for a 
more conservative choice of a 8 s in connection with zeros * 

The solution obviously meets static specifications 
because of the introduction of the feed forward integration 
compensation* 

8*5 Summary of the Application of the Technique to this Problem 
The solution of the problem proceeded along straight- 
forward lines until it was discovered that a Position Control 
System could not provide acceptable following error* It was 
then necessary to stipulate compensation which would give the 
system Velocity Control characteristics* With this compen- 
sation inserted, the solution was then recommenced* The 
necessity of another compensator became apparent and its 
characteristics were determined* 

It is interesting to note the similarity between the 
compensation determined by this solution and that provided 
by Chestnut and Mayer Both involve feedback of first 

and second derivatives of the output* However, this solution 

2 

replaces their elaborate filter, fts ^s+a) , with 

(s+b) (s+c) (s+d) 

feed forward through the function K and feedback of a signal 

s 

proportional to the output quantity* 
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CHAPTER 9 

CONCLUSIONS AND RECOMMENDATIONS 

9.1 Conclusions . Application of step-by-step procedures 
developed in this thesis resulted in successfully designing 
compensation for a system which, though not overly complex, 
was subject to somewhat severe specif ications . Alternate 
methods of compensation design may have resulted in equally 
simple comepnsation but, by their trial-and -error nature, 
must be assumed to be less straightforward <> It it there- 
fore concluded that the proposed compensation technique is 
valid and presents the designer with a powerful and labor- 
saving tool. 

9*2 Recommendations <> Several areas for further study present 
themselves upon consideration of the results of this thesis 
and of the steps by which these results were obtained. 

Specific areas requiring further work are: 

a. More exact determination of the deviation from 
second-order response caused by the existence of 
zeros and additional roots. 

b. Careful consideration of the step-by-step procedure 
to eliminate possible duplication of effort when order 
raising compensation is used. 

Areas of research which should provide beneficial results ares 
a. Examination of the means by which specifications can 
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be related to dominant third-order response* 

b. Consideration of creating non-dominant oscillatory- 
modes at harmonics of the dominant mode, in order to 
better reproduce the wave form of the input. 

c. Examination of the means by which this method may 
be extended to include treatment of other than posi- 
tional servomechanisms. 
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APPENDIX A 



A DETAILED EXAMINATION OF THE EFFECT OF a ON 
THE DEVIATION OF A THIRD ORDER RESPONSE FROM THAT 
OF A SECOND ORDER 



To examine the effect of on the deviation of the third 
order response from the second order, consider the variation 
of the two systems if a unit step is applied to each. 



Second Order: Out(s) = p 5— 

s(s +2Co> s+co ) 
n n 



(Ad) 



n 



where Phase lag of ideal (2nd order) response = 
?2 



- arctan 



< 



Third Order: Out(s) = 



a l m n 



s ( s+ a, <o ) ( s 2 +2 Ca>_ s+a) 2 ) 
in n n 



(A.3) 



where a^co n is included to make steady state values equal. 



-I -aoo t 

Out(t) = — p p e 9 — + 

oo n co n ^[a^-2Cd+l] 

Q 1 e" C “n t sinG/l-< 2 (o n t+l]! a 
co n 2 [ ( 1 - C 2 ) ( a 2 - 2?: a + 1 )] 1 / 2 



(A.*+) 
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where lf a = Phase lag of actual (3rd order) response = 

, 2 / 2 

- arctan vl-C - arctan vl-K . 

-C a-C 

Define deviation as follows: 

D = actual - ideal (A. 5) 

ideal 

Substitution of^equations (A. 2) and (A. 4) into equation 
(A. 5) yields: 



e - ( a -O (Oflt + a sln(v^C 2 co n t^ a ) sln( ^2 



D = 



[a -2Ca+l] 



[a -2ta+}] 



C0 n t+ ^i’ ) 



/ 1 < 2 + sinC^T^ -2 co n t + TjT ;L ) 



(A. 6) 



This deviation is a function of t,a, and C of the dominant 
roots. Ideally we would like to find the maximum deviation, 
no matter when it occurs. This, however, is a great under- 
taking which is beyond the scope of this paper. To simplify 
the problem we have picked the time corresponding to the peak 
overshoot of the second order response: 



T = 
P 



~2 + arctan 



/i< 2 



Substituting this time into equation (A. 6): 



/r-v 2 

1-C e 



-(a-OooT 
n p 



a sin(^+arctan^f|- 



D - f q -2Cq+1 1 






r a 2 ~2Ca+ll 1</2 



-1 



n P + 1 



(A. 7) 
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but Sin (g+arctan ') = 



[a 2 -2Ca+l) 1//2 



so: 



D = 




n P + a£+l 



(A. 8) 



[lVi-C* e" J “n T P][a 2 -2Ca+l] 



This equation is plotted in figure A.l. Figure A .2 represents 
the result of similar consideration of the effect of a zero 
at ao> n on a second order response. 

Examination of these curves will yield an appropriate 
value of a to insure the desired degree of dominance of the 
second order response. 
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APPENDIX B 



THE MODEL CHARACTERISTIC EQUATIONS 



B-l Second-Order Equations 

2 

s + a n s + = 0 

1 o 



Minimum Model: 



a l = 2C<0 n 

a = <o 2 
o n 



Working Model: 



= a> 
o n 



B-2 Third-Order Equations 
2 

s + a 0 s + a n s + = 0 



2^ ' a l 

Minimum Model 



a 2 = 

a 1 = (l+2Ca 1 )o> n 2 



a o = (0 1 } V 



a 2 = (2C+a^)o> n 

a, = (1+2CCU )« 2 
1 in 



a o = (a l )c V 



Working Model: 
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a 3 " 


2 ^n 


II 

CM 

cd 


(l+a 1 2 )a> n 


II 

1— 1 
cd 


< 2 <°i 2) % 


II 

o 

cd 


, 2. 4 

(o l )co n 



B-3 Fourth-Order Equations 
L ^ 2 

s + a^s^ + a 2 s + a-^s + a Q = 0 
Minimum Model: 



Working Model (Non-dominant roots complex) 
a 3 = (2^+2^)^ 

a 2 = (l +l+a i^i^2 +a l 2 ^ co n 2 

a l = V'lVXl'lK 3 

a o = (a l 2)co n 4 

Working Model (Non-dominant roots real) 
a^ = (2C+P 1 )co n 

a 2 = (1+2CP 1 +P 2 )£0 n 2 



a ± = (P 1 +2CP 2 )a) n J 
a o = (P 2 )to n 4 



where 



P 1 = Q l +0 2 



>2 = “l°2 



B-4 Fifth-Order Equations 



r* + a^s^ + a ^ 
Minimum Model: 



s + a^s + a 3 s + a 2 s + a-^s + a Q = 0 



a 4 = (2!: i +a 2 )co n 
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B-b Fifth Order Equations (cont'd) 

Minimum Model (cont'd) 

a 3 = (l+o 1 2 +2C 1 a 2 )w n 2 
a 2 = (2t; i a 1 2 + a 2+ a 1 2 a2) ( o n 3 

a l = («i 2 +2Ci-i 2 a 2 )co n ‘ + 

a Q = (a! 2 ^)^ 5 

Working Model (One pair non-dominant roots complex) 

a 4 = (^i + ^ a i^2 +a 2 ^n 

a 3 = (l+4a 1 !: i C 2 +a 1 2 +2!: i a 2 +2a 1 a 2 C 2 )a3 n ‘ i 

a 2 = (2a 1 C 2 +2C 1 a 1 2 +a 2 +4a 1 o 2 C 1 !: 2 +a 1 a 2 )co n 3 
2 2 If 

a l = (a-L +2a 1 a 2 C 2 +2<x 1 

a - (a.^a 0 )co ^ 
o v 1 2 n 

Working Model (Non -dominant roots real) 

a 4 = (2!: 1 +P l )<0 n 
a 3 = (l+2C 1 P 1 +P 2 )co n 2 



a 2 = ^1 +2 ^ P 2 +P 2 )£ “n 3 
a l = (P 2 +2!: i P 3 )a) n 1+ 
a o = (P 3 )co n 5 



P 1 = °1 + °2 + “3 

p 2 = ° 1 a 2 + a 2 a 3 + a 3 a l 

P 3 = a l ° 2 a 3 



where 
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B-5 Sixth-Order Equations ; 

+ a^s^ + + a 3 s ^ + a 2 s ^ + a i s + a ~ = 0 

Minimum Model: 



a 5 " (2! *l )a n 



= (l+a 1 2 +a 2 2 )o5 n 2 

a 3 = ( 2 ^l a l 2+2 ^l a 2 2 ^n 3 
, 2^ 2^ 2 2> 4 

a 2 (Q 1 + °2 + °1 a 2 )a) n 

2 2 , 5 

a l = (2C l a l a 2 )<0 n 
, 2 2 . 6 
a o = (a l a 2 )£0 n 



Working Model (Both pairs non-dominant roots complex) 

a^ = (2C^ + 2a^C2 + 2tt2^3^ a) n 

a 4 = (1+4a l !: 1 C 2 +a l 2+4!: 1 ' : 3 a 2 +4!: 2 !: 3 a 1 “ 2 +a 2 2 )“ n 2 
a 3 = (2a 1 C 2 +2t: i a 1 2 +2a 2 !: 3 +8a 1 a 2 !: i t: 2 !: 3 +2!: 3 a 1 2 » 2 

+2 h a 2 2+2a 3. a 2 2t: 2 )M n 

a 2 = ^“l +lfo i a 2 , ’2 ! *3 +l<a l a 2 ! "l ! ’ 3 +a 2 +1+, 'l ! ’2 a l a 2 

a l = ( 2 ^3 a l 2a 2 +2a l a 2^2 +2a l 2o 2 2 ^l^n 5 
, 2 2 . 6 
a o = ( °1 °2 )a >n 

Working Model (One pair non-dominant roots real) 
a ^ = (2C 1 +2a 1 C 2 +P 1 )co n 

= (l+>+a 1 C 1 !; 2 +a 1 2 +2!: i P 1 +2!: 2 a 1 P 1 +P 2 )oo n 2 
a 3 = (2a 1 C 2+ 2!: i a 1 2 + P 1 ^a 1 P l! : ;L !: 2+ a 1 2 p i+ 2 ! : i P 2 



+2a l P 2 !: 2 ) “n 
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B-5 Sixth-Order Equations (cont'd) : 



Working Model (One pair non-dominant roots real) 

i- a. 

5 



a 2 = (« 1 2 +2a i p i! : 2+ 2 ! ; i a 1 2 p i+ p 2+ 4a i p 2! : i! : 2+ a 1 2 p 2 )c0 n ‘ + 



a l = (a l 2p l +2 “l P 2 ! -2 +2a l 2p 2 t *l^“n y 
a o = (a l 2p 2 )to n 6 



where 



3 1 = a 2 +a 3 



*2 = a 2“3 



Working Model (All non-dominant roots real) 

a 5 = (2C 1 +P l )(0 n 
a 4 = (1+2!: 1 P 1 +P 2 )» n 2 

a 3 = (P 1 + 2C 1 P 2+ P 3 )co n 3 

a 2 = ^2 +2 ^3 + ^ )l °n 
a l = ^ 3 +2C l P V )a3 n 5 
a o = (P 4 )£ “n 6 

where 

P 1 = Q i +a 2 +a 3 +a 4 

P 2 = a 1 a 2 + a 1 a 3 + a 1 a 1 ++ a 2 a 3 + a 2 a lf - t - a 3 a 4 

P 3 = flV3 +a l a 2 a 4 +a 2 a 3V a l a 3\ 

p 4 = a 1 a 2 a 3 a l^ 



B-6 Seventh-Order Equations ? 

s ^ + a^s 6 + a^s^ + a^s^ + a ^s^ + ^s 2 + a n s + a^ = 



0 



9 ^ 

Minimum Model; 



a 6 = 


(2C l +a 3 )<0 n 


a 5 = 


(l+'a 1 2 +a 2 2 +2t: i o 2 )to n 2 


a 4 = 


(2C 1 a 1 2 +2C 1 a 2 2 +a3+a 1 2a 3 + a 2 2a 3 )a, n 3 


a 3 = 


(«i 2 + . 2 2 + «i 2 . 2 2 +2 c 1 -i 2 . 3 + 2 «:i« 2 2 .3 ) ^ 4 


a 2 = 


(2!: i a l 2a 2 2+a l 2Q 3 +a 2 2o 3 +a l 2a 2 2a 3 )a) n 5 


a l = 




a o 





Working Model (Some roots complex) 



a 6 = 


(A^+a)co n 


a 5 = 


(A 1+ +aA^)03 n 


a l+ = 


(A^+aA^)a) n 


a 3 = 


(A 2 +aA^)a) n 


a 2 = 


(A 1 +aA 2 )(O n 


a l = 


( A 0 +aA 1 )“ n 


a o = 


(aA o )co n 


where A refers to 
n 


coefficients of same type model in 


sixth order eqns. 




Working Model (All 


non-dominant roots real) 


a 6 = 


(2 h +P l )w n 


a 5 = 


(1+2C 1 P 1 +P 2 )oo n 2 


a 4 = 


(P l +2 h P 2 + P 3 )( “n 3 


a 3 = 


(p 2 +2C l P 3 +p 4- )<0 n lf 
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a 2 = (P 3 +2C 1^ +P 5 )w n 5 

a i = 

a o = (P 5 )co n 7 

where 

p i = a i +a 2 +0 3 + V a 5 

p 2 = a i a 2 +Q i a 3 +a lV Q i a 5 +a 2 a 3 +Q 2V a 2°5 +Q 3V a 3 a 5 + % a 5 
P 3 = a 1 . 2 a 3 + « 1 « 2 V a 1 « 2 . 5 + « 2 « 3 a i ; + « 2 a 3 a 5 + c 3 a lf a 5 

+a l a 3V a i a 3 Q 5 +a i\ a 5 +a 2 \ a 5 
P 4 = a l a 2 a 3% +a l a 2 a 3°5 +a 2 a 3\ a 5 +? l a 3 a 4 a 5 +a l a 2\ a 5 

P 5 = o 1 a 2 a 3 a l+ a 5 
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